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We present calculations of the interplane charge dynamics in the normal state of cuprate super- 
conductors within the valence-bond dynamical mean-field theory. We show that by varying the 
hole doping, the c-axis optical conductivity and resistivity dramatically change character, going 
from metallic-like at large doping to insulating- like at low-doping. We establish a clear connection 
between the behavior of the c-axis optical and transport properties and the destruction of coherent 
quasiparticles as the pseudogap opens in the antinodal region of the Brillouin zone at low doping. 
We show that our results are in good agreement with spectroscopic and optical experiments. 

PACS numbers: 71.27.-|~a, 74.72.-h, 74.25.-q 



I. INTRODUCTION 



Even after many years of investigation, cuprate super- 
conductors continue to be at the center of intense ex- 
perimental and theoretical interest. A key phenomenon 
which needs to be understood to reveal the nature of su- 
perconductivity in the cuprates is the onset of strong 
momentum-space differentiation in the (hole-) under- 
doped normal state. This region of the phase diagram is 
characterized by the suppression of quasiparticle excita- 
tions in the antinodal region of the Brillouin zone and the 
opening of a pseudogap yielding Fermi circs, 0jS observed, 
e.g. in angle-resolved photoemission (ARPES) experi- 
mentsii Other distinctive properties of the cuprates are 
seen in the extreme anisotropy of the charge dynamics.—^— 
The in-plane {pab) and interplane resistivity (pc) have dif- 
ferent temperature dependence. While pab decreases with 
decreasing temperatures at all doping levels, pc is sensi- 
tive to the doping. Its behavior is essentially insulating 
at low doping, semi-metallic at intermediate doping and 
metallic in the overdoped regime. The peculiar proper- 
ties of the interplane charge dynamics have also been ob- 
served in optical studiesi^"— They reveal that the c-axis 
optical conductivity Oc{^) shows no sharp Drude peak 
at low frequencies in the undcrdopcd normal state. In- 
stead, i7c(ri) has a rather flat frequency-dependence at 
temperatures above the pseudogap temperature T* . As 
the temperature is reduced below T* , low-energy spectral 
weight is transferred to higher energies and crc(O) dis- 
plays a gaplike depression as 0. Only in the highly 
overdoped region does (Tc(f^) show evidence of emerging 
coherence 4 

In this article, we address the interplane charge dynam- 
ics of cuprate superconductors within the valence-bond 
dynamical mean-field theory (VB-DMFT) introduced in 
Refs. [13 and [3. Our results capture many salient fea- 
tures found in experiments for the c-axis resistivity and 



optical conductivity. Specifically, we show that the open- 
ing of a pseudogap in the antinodal region explains the 
incoherent behavior of the interplane transport. 




FIG. 1. (Color online) The two patches dividing the Brillouin 
zone. The line shows a non-interacting Fermi surface for the 
dispersion of Eq. [1] The central (red) patch covers the 
nodal region of the Fermi surface, while the border (blue) 
patch covers the antinodal region. 



The VB-DMFT is a minimal cluster extension of the 
dynamical mean-field theoryi^ji^ that aims at describ- 
ing momentum-space differentiation together with Mott 
physics. We apply it to the Hubbard model on a square 
lattice defined by the Hamiltonian 



(1) 



k,cr 



Ek = —2t[cos{kx) + cos(ky)] — 4:t' cos(kj,) cos(ky), 

where t and t' are the nearest and next-nearest neighbor 
hopping and U is the onsite Coulomb repulsion. In the 
following, we use U/t = 10 and t' /t = —0.3, which are 
values commonly used for modeling hole-doped cuprates 
in a single-band framework. This value of U is larger than 
the critical value for the Mott transition in the undoped 
model, so that we deal with a doped Mott insulator. 
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FIG. 2. (Color online) Left panel: The c-axis optical conductivity ac{Q) calculated within VB-DMFT for three doping levels. 
(Jc is displayed in units of a* as defined in the text {a* is of order 50 Q~^cm~^ for YBa2Cu30y). Frequency is normalized to 
the half-bandwidth D ~ leV — 8000cm~^. Rightpanel: Experimental data for the c-axis optical conductivity of YBa2Cu30j,. 
The data for YBa2Cu307.oo is taken from Ref. where the phonon contribution was subtracted by fitting to 5 Lorentzian 
oscillators. The data for YBa2Cu3 06.95 and YBa2Cu30c.67 are taken from Refs. [iB andllTl. 



VB-DMFT is based on an effective two-impurity model 
embedded in a self-consistent bath. Within this simple 
approach involving only two degrees of freedom, the lat- 
tice self-energy is approximated to be piecewise constant 
over two patches of the Brillouin zone. These two patches 
are shown in Fig. [TJ Inspired by the phenomenology of 
cuprate superconductors, their shape is chosen in such 
a way that the central (red) patch covers the nodal re- 
gion of the Fermi surface, while the border (blue) patch 
covers the antinodal region. With this prescription, one 
degree of freedom of the underlying two-impurity model 
is associated with the physics of nodal quasiparticles and 
the other degree of freedom with antinodal quasiparti- 
cles. Thank to the relative simplicity of the VB-DMFT 
approach, it is possible to make very efficient calculations 
using continuous-time quantum Monte Carlo and obtain 
accurate spectral functions A(k,uj). More details con- 
cerning this procedure are given in Appendix |^ and in 
Refs. [ij and tli 

In the VB-DMFT, the formation of Fermi arcs is de- 
scribed by a selective metal-insulator transition in mo- 
mentum spaceM^ Below a doping ^ 16%, the degree 
of freedom describing the antinodal regions becomes in- 
sulating, while that associated to the nodal quasiparti- 
cles remains metallic. The orbital-selective mechanism 
responsible for the pseudogap has also been confirmed 



in studies involving larger clustersi^"— In Refs. [13 and 
[Tsl . we used the VB-DMFT to compute tunneling and 
ARPES spectra in good agreement with experiments. 



II. INTERPLANE OPTICAL CONDUCTIVITY 

We first compute the frequency-dependent c-axis opti- 
cal conductivity ad^), given by 

nab / s I 

(2) 

l5]4(k)A(k,c.)A(k,^] + c.), 

k 

where / is the Fermi function, A(k, uj) the in-plane spec- 
tral function, the number of lattice sites, e the elec- 
tronic charge, a, b the in-plane lattice constants, c the 
interplane distance and t±(k) = to(cos(k2;) — cos(kj,))^ is 
the interplane tunneling matrix element)^— Note that 
in this expression t±{'k) has a strong k-dependence, with 
contributions stemming mainly from the antinodal re- 
gion of the Brillouin zone. For convenience, we will 
express energies in units of the half-bandwidth D of 
the electronic dispersion, and the optical conductivity 
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in units of a* = 2e'^ct'^/habD^ . In YBa2Cu30y com- 
pounds, D ^ leV ^ SOOOcni^^ and a* is of order 
cr* - SOrj-^cni-l. 

In the left panel of Fig. [31 we display the computed 
crc(fi) for three levels of hole-doping and several tem- 
peratures. Our results show three distinctive behaviors. 
At high doping S > 16%, the conductivity displays a 
metallic-like behavior, with the build-up of a Drude-like 
peak as the temperature is decreased. Note that as the 
peak increases additional spectral weight appears at low- 
energy. At low doping 5 < 10%, (Jci^) is characterized 
by a gaplikc depression at low frequencies where spectral 
weight is suppressed with decreasing temperature. The 
width of the gap when it opens at high temperature is 
^ 0.15-D and remains approximately the same as the tem- 
perature is lowered. Note that the spectral weight that is 
lost in the gap is redistributed over a wide range of ener- 
gies. The appearance of the depression in the spectra can 
be directly linked to the formation of a pseudogap in the 
antinodal regionJ^ii^ Indeed, the matrix element t± ap- 
pearing in the expression of the optical conductivity ^ 
essentially probes the regionS^ close to (±7r, 0), (0, ±7r) so 
that a loss of coherent antinodal quasiparticles results in 
a loss of low-energy spectral weight in the c-axis opti- 
cal conductivity. In Refs. [13 and Tisl it has been shown 
that in a zero-temperature analysis of VB-DMFT, co- 
herent quasiparticles disappear in the antinodal region 
at a doping ~ 16%. This is consistent with CTc showing 
a depression only for doping levels below ~ 16%. At 
intermediate doping, ad^) has a mixed character (see 
left panel of Fig. [2] for (5 ~ 12%). Starting from high 
temperatures, the low-energy conductivity first displays 
metallic behavior with an increase in CTc as the temper- 
ature is lowered. At temperatures T < D/50, however, 
a gap starts to appear and low-energy spectral weight is 
suppressed. 

The right panel of Fig. [2] displays the experimental 
interplane optical conductivity of YBa2Cu3 0y obtained 
in Refs. d, [l^, and [l3- At large doping levels, the data 
shows a build-up of a Drude peak consistent with our 
theoretical calculation. One must note however that the 
interplane optical conductivity of YBa2Cu307 displays 
a crossing between the spectra at different temperatures 
and that the spectral weight transfer is not as clear as 
in the theoretical curves (left panel) for S = 16%. In 
fact, only at doping levels larger than S ^ 25% do the 
computed ac eventually behave more metallic-like, with 
a sharp Drude peak at low energies. This behavior is 
also consistent with other experiments, e.g. on overdopcd 
YBa2Cu30y, sec Fig. 1 of Ref. [2^. At low doping, the 
spectra of YBa2Cu306.67 display the opening of a gaplikc 
depression, which is quite well captured by our theoret- 
ical results. Remarkably, CTc is typically a factor of 10 
smaller that in YBa2Cu307, a factor quantitatively sim- 
ilar to what we obtain theoretically by comparing our 
results for 6 — 16% and 5 = 8%. At intermediate doping 
between these two limits, YBa2Cu3 06.95 displays spectra 
which depend on temperature in a non-monotonous man- 



ner, first increasing upon cooling and then decreasing at 
lower temperatures. This behavior is indeed observed in 
our theoretical results at 12% doping. All these observa- 
tions show that our calculations agree qualitatively, and 
to a certain extent quantitatively, with optical conductiv- 
ity measurements on underdoped cuprates (see also e.g. 
Ref. i. Fig. 2 in Ref. or Fig. 1 in Ref. [13). 



III. INTERPLANE RESISTIVITY 

From the £7 = extrapolation of the optical conductiv- 
ity, we obtain the c-axis resistivity pc = l/ad^ = 0) dis- 
played in Fig.|3]for different doping levels, as a function of 
temperature, pc is displayed in units of — 1/a* which 
takes the value ^ 0.02f2cm in YBa2Cu3 0j,. As already 
discussed for the optical conductivity, the properties of 
Pc strongly depend on the doping level. At low doping 
{S = 8%), Pc{T) has an insulating behavior with a large 
increase as the temperature is lowered. This increase is 
a direct consequence of the opening of the pseudogap 
in the antinodal region of the Brillouin zone. As the 
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FIG. 3. (Color online) The c-axis resistivity pc calculated 
within the VB-DMFT for four doping levels as a function of 
temperature. 

doping level is increased (see 6 = 12%), pc{T) starts to 
display a mixed, non monotonous behavior: it has a min- 
imum as it changes from metallic-like to insulating-like 
as T decreases. This behavior is indeed observed in inter- 
plane resistivity measurements on La2_i^Sr^Cu04 (Fig. 3 
in Ref.i, Fig. 2 in Ref. [3), YBa2Cu3 0y (Fig. 1 in Refjl 
or Bi2Sr2CaCu20j, thin films (inset of Fig. 3 in Ref. . 
As we go to higher doping, Pc{T) eventually acquires a 
metallic-like behavior with dpc/dT > 0. Note that this 
metallic behavior emerges from the transfer of spectral 
weight to lower energies as can be seen from Fig.[51 More- 
over, even at the largest doping level S = 25%, Pc{T) 
does not recover the standard Fermi-liquid behavior with 
Pc{T) ^ (closer inspection instead shows that Pc{T) 
roughly behaves as T" with a k, 1.6). One might need 
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to consider even larger doping levels to recover a be- 
havior as seen, e.g. in heavily doped La2_2;Sr2;Cu04i^ 



IV. CONCLUSIONS 

To conclude, in this article we have computed the inter- 
plane charge dynamics within the VB-DMFT framework 
in order to address the properties of the normal-state of 
high- Tc copper-oxide superconductors. Our results show 
that the intcrplane charge dynamics of the normal state is 
characterized by three regimes, (i) An overdoped regime 
where both the c-axis optical conductivity and resistivity 
show a metallic-like behavior, (ii) An underdoped regime 
where the destruction of coherent quasiparticles in the 
antinodal region induces a gaplike depression in CTc and 
a resulting insulating-like interlayer resistivity which in- 
creases as temperature is lowered, (iii) An intermediate- 
doping mixed regime with first the buildup of a coher- 
ence peak in CTc as the temperature is decreased, followed 
by a loss of low-energy spectral weight as the pseudogap 
opens at lower temperatures. As a consequence, pc has a 
non-monotonic behavior with a minimum at intermedi- 
ate temperature. Our calculations are in good agreement 
with spectroscopic and transport experiments. 



ACKNOWLEDGMENTS 

We are grateful to Hclcnc Raffy for discussions about 
her experimental results on c-axis resistivity. We ac- 
knowledge the support of the Partner University Fund 
(PUF-FACE), the Ecole Polytechniquc-EDF chair on 
Sustainable Energies, the National Science Foundation 
(NSF-Materials World Network program and NSF-DMR- 
0906943), and the Agence Nationale de la Recherche 
(ANR grant ECCE). This work was performed us- 
ing HPC resources from GENCI-CCRT (grant 2009- 
t2009056112). We also acknowledge the Kavh Institute 
for Theoretical Physics where part of this work was com- 
pleted (under grant NSF-PHY05-51164). DNB is sup- 
ported by the NSF. 

Note added: After initial submission of this 
manuscript, Lin, Gull and Millis (arXiv: 1004,2999) re- 
ported calculations of c-axis optical conductivity in 8-site 
cluster DMFT. When comparison is possible, consistency 
between these results and our VB-DMFT (2-site) results 
is observed. 



DMFT, the effective model that one has to solve is a 
two-impurity Anderson mode l^^i^^ . The even (-|-) combi- 
nation of the two impurity orbitals is associated with the 
nodal part of the Brillouin zone, while the odd (-) com- 
bination is associated to the antinode. The Anderson 
impurity model is solved using a continuous-time quan- 
tum Monte Carlo algorithm (CTQMC)^i2a that yields 
the even (rcsp. odd) orbital self-energy E_|_(iw„) (resp. 
I]_ {iujn)) at the Matsubara frequencies At this stage 
one can readily get the quasiparticle residues by carefully 
extrapolating Ti±(iijjn) at zero frequency 



Z±= 1 



dl'niT,±{iLOn) 



duj„ 



(Al) 



The inverse quasiparticle lifetime is then given by 



r± = -Z±llJlT,±{iuJn) 



(A2) 



Finally, the nodal (Fn) and antinodal (Fan) inverse life- 
times are associated with r+ and F_ respectively. In 
order to compute the optical conductivity and the re- 
sistivities it is necessary to have access to the in-plane 
spectral function A(k, w). A first step is to analytically 
continue the self-energies Y,± on the real-frequency axis, 
which is achieved using Pade approximantsi^ This is 
made possible thank to the extremely accurate, low-noise 
imaginary-time data. A second step is to interpolate the 
self-energy to obtain a k-dcpendent I](k, w) from which 
the spectral function will be computed. Refs. [13 and 
fill have shown that an efficient interpolation procedure 
is the Af -interpolation.— 1^ The self-energy is obtained 
with S(k, Lu) = Lu + p, — M(k, w)~^, where the cumulant 
M(k, w) is given by 



M(k,a;) 



a+(k) 



uj + fi — w + /.j — S_ (oj) 



(A3) 



with Q!±(k) = i{l ± i[cos(ka;) -|- cos(kj,)]}. Finally, the 
spectral function is computed with 



Appendix A: Details about the VB-DMFT approach 

Let us give here some technical details and briefiy de- 
scribe how our results were obtained. Within the VB- 



A{k,uj) 



w + /i - Ek - S(k, w) 



where is the square lattice dispersion. 



(A4) 
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